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After a brief presentation of the MACRO detector we discuss the updated data on 
atmospheric muon neutrinos, and the interpretation in terms of neutrino oscilla- 
tions. 



1 Introduction 

MACRO is a multipurpose underground detector designed to search for rare 
events in the cosmic radiation. It performs measurements in areas of astro- 
physics, nuclear, particle and cosmic ray physics. In this paper we shall discuss 
the measurement of the atmospheric muon neutrino flux in the energy region 
from a few GeV up to a few TeV, and the interpretation of the MACRO re- 
sults in terms of neutrino oscillations. The detector has global dimensions of 
12 x 9.3 x 76.5 m 3 and provides a total acceptance to an isotropic flux of par- 
ticles of ~ 10, 000 m 2 sr. The. total mass is ~ 5300 t. A cross section of the 
detector is shown in Fig. [IJaLy. It has three sub-detectors: liquid scintillation 
counters, limited streamer tubes and nuclear track detectors. The mean rock 
depth of the overburden is ~ 3700 m.w.e. The average residual energy and the 
muon flux at the lab depth are <~ 310 GeV and ~ 1 rcT 2 h _1 , respectively. 

2 MACRO as atmospheric v„ detector 

Primary cosmic rays produce in the upper atmosphere pions and kaons, which 
decay, ir — ► [iv, K — > [iv\ and fj, — * eVf.v^. Fig. |l]a shows the three different 
topologies of neutrino events detected by MACRO: up throughgoing muonsH, 
semicontained upgoing muons Q.U) and up stopping muons + semicontained 
downgoing muons (UGS + ID) a. The parent energy spectra for the three 
event topologies were computed with Monte Carlo (MC) methods, and pre- 
sented in Fig. ^p. The number of events measured and expected for the three 
topologies is given in Table |l|. Sources of background and systematic effects 
were studied in detail in H Q and found to be negligible. The neutrino energies 
involved in all event topologies are such that the /i„ flux is up-down symmetric. 

High energy sample. The up throughgoing muons come from v^s in- 
teracting in the rock below the detector; they are identified using the streamer 
tube system (for tracking) and the scintillator system (for time-of-flight (ToF) 

a Invited talk at the Third International Workshop on New Worlds in Astroparticle Physics 
1-3 September 2000, University of the Algarve. Faro, Portugal 



1 




Figure 1: a) Event topologies induced by interactions in and around MACRO. The 
black boxes are scintillator hits. The track direction and versus are measured by the 
streamer tubes and time of flight, b) Neutrino energy distribution for the three event 
topologies detected by MACRO 
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Table 1: Event summary for the MACRO atmospheric neutrino flux analyses. The 
data correspond to a lifetime of ~ 5 y with the full detector, after background 
corrections [2,3]. The ratios R = Data/MC are relative to MC expectations assuming 
no oscillations. 

measurement). A rejection factor of at least 10 5 is needed in order to separate 
the up- going muons from the large background of down-going muons □. The 
median energy is E v ~ 50 GeV , and muons with > 1 GeV cross the 
whole detector. Fig. ||a shows the zenith angle distribution of the measured 
flux of lip throughgoing muons. The data have been compared with MC simu- 
lations! 2 ]; the expected angular distributions for no oscillations is shown in Fig. 

as a dashed line, and for a — + v T oscillated flux with sin 2 26* = 1 and 
Am 2 = 0.0025 eV 2 , as a solid line. The total theoretical uncertainty on the 
expected muon flux, adding in quadrature the errors from neutrino flux, cross 
section and muon propagation, is 17 %; it is mainly a scale error that does 
not change the shape of the angular distribution. The ratio of the observed 
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Figure 2: (a) Measured flux (points) of the up throughgoing muons plotted vs. zenith 
angle Q. The solid line is the fit to an oscillated muon flux, obtaining maximum 
mixing and Am 2 = 0.0025 eV 2 . (b) Confidence regions for —* v T oscillations at 
the 90 % CL, obtained from MACRO high energy [2] and low energy [3] data. 



number of events to the expectation without oscillations is given in Table |. 

Low energy samples. The upgoing semicontained muons (IU, in Fig. 
|l|a) come from interactions inside the lower apparatus. Since two scintilla- 
tion counters are intercepted, the ToF is applied to identify the upward going 
muons. The average parent neutrino energy for these events is 4.2 GeV. The 
upgoing stopping muons (UGS) are due to external interactions yielding 
upgoing muon tracks stopping in the detector; the semicontained downgoing 
muons (ID) are due to induced downgoing tracks with vertices in the lower 
MACRO. The events were selected by means of topological criteria; the lack 
of time information prevents to distinguish the two sub samples, for which an 
almost equal number of events is expected. The average neutrino energy for 
these events is ~ 3.5 GeV . The number of events and the angular distributions 
are compared with the predictionsB in Fig. || and Table |l|. 

3 Interpretation in terms of neutrino oscillations 

We interpreted the reduction on the detected number of events (Table 1) and 
the deformation of the zenith angle distribution of the up-throughgoing events 
(Fig. 2a) as a consequence of disappearance. 

The weak flavor eigenstates (u e , , v T ) are relevant for the production of 
atmospheric neutrinos, via 7r — > (Xfu, K— > [iv^ and [i — > ev^^ decays. For mas- 
sive neutrinos one has to consider the mass eigenstates (i>\, 1/2, v z) which are 
relevant for propagation. The flavor eigenstates are linear combinations of the 
mass eigenstates, vi = ^ m=1 Ui m v m . In the simplest two flavors oscillations 
scenario, only one mixing angle 9 is involved, and the rate of disappearance 
is given by 

Piy^ -> v^) ~ 1 - sm 2 20 sin 2 (1.27 Am 2 L/E v ) (1) 
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Figure 3: Measured (points) and expected number of low energy muon neutrino events 
versus zenith angle. IU: up semicontained. ID+UGS: up stopping plus down semi- 
contained muons. The shadowed regions are the predictions without oscillations; the 
dotted lines are the predictions assuming neutrino oscillations with the parameters 
obtained from the up throughgoing sample. 



Am 2 = m 2 „ 3 — m 2 V2 , L — distance (in km) from the decay point to the 
interaction point, E u is the neutrino energy, in GeV. 

This simple relation could be modified .when neutrinos propagate with a 
flavor-dependent interaction through mattertl This case is important, because 
the matter effect could help to discriminate between oscillation in different 
neutrino channels. The matter effect in the Earth could enhance (in a range 
of values of the oscillations parameters) the disappearance effect through 
"resonances" (MSW effect). This happens only for — > v e and for — > v s 
channels [y s = sterile neutrino), while for — » v T there is no matter effect 
{Vfx and v T have the same weak potential). 

In the scenario described by eq. 1, relatively fewer up throughgoing 
muons are expected near the vertical (cosQ = —1) than near the horizon- 
tal (cosO — 0), due to the longer path length of neutrinos from production to 
observation, in a wide range of Am 2 values. The probabilities to obtain the 
number of events in Table 1 and the angular distribution observed in Fig. 2a 
have been calculated vs. oscillation parameter values. For — > v r oscillations 
the maximum probability is 57 %, for Am 2 = 0.0025 eV 2 , sin 2 29 = 1 (com- 
bination from the shape of the angular distribution and the reduction in the 
number of events). Fig. ||b shows the 90 % CL regions for — ► v T . 

The probability for no-oscillations is 0.4%. The maximum probability for 

— > v s oscillations is 15%. Another way to discriminate between the 
oscillation in v T or in v s is to study the angular distribution in two bins. The 
statistical significance is higher □ with respect to the 10-bins fit of the angular 
distribution. From this ratio, the — > v s hypothesis is disfavored at level of 
2.5cr compared to — > v^a. 
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From our results on the high energy sample, a ~ 50% reduction in the 
flux of the up stopping events (UGS) as well of semicontaincd upgoing muons 
(IU) is expected, with no distortion in the shape of the angular distribution. 
No reduction is instead expected for the semicontained downgoing events (ID, 
coming from neutrinos with path lengths of ~ 15 km). This is what is observed 
m Fig. |. The allowed region in the Am 2 , sin 2 29 plane is shown in Fig. ^b. 

4 Up throughgoing muon energy separation 

The up throughgoing muon sample is induced by atmospheric in a wide 
energy range (for 90% of events 5< E v < 700 GeV, and 2< <300 GeV, see 
Fig. lb). If the observed event reduction is due to neutrino oscillations, the 
first minimum of eq. 1 occurs at E u ~ 10 — 30 GeV (for neutrinos from the 
nadir direction and oscillation parameters in the rang c of Fig. |b). Thus, the 
f M disappearance mechanism works on (relatively) lower energy neutrinos. As 
a consequence, we expect (from MC) a ~ 50% reduction of the up through- 
going muons with E^^IO GeV , and only a ~ 20% reduction for muons with 
£ M >10 GeV. 

Muons traversing MACRO are deflected by many small-angle scatters, 
the bulk of which are Coulomb scatterings from the nuclei and the atomic 
electrons (multiple Coulomb scattering). As an effect, relatively low energy 
up throughgoing muons are (slightly) deviated from a straight-line direction 
inside the detector. We evaluated that this deviation is of the order of few 
centimeters, up to E^ ~ 10 GeV, and it exceed the intrinsic spatial resolution 
of the MACRO streamer tube cells. 

A selection cut, based on the deviation of the muon sagitta from the 
straight line when the muon cross the whole detector, was defined. Each 
up throughgoing muon is classified as: Low= muon with a deviation from the 
straight line (from MC, lower energy muons with median E^ ow = 11 GeV); 
High= muon with no measurable deviation (muons with median energy E^ 9 — 
52 GeV). The remaining events with ambiguous deviation, and which does not 
satisfy the high or the low energy cuts, were classified as Medium. From MC, 
they have an intermediate energy Eff mm = 33 GeV . The same cuts were 
then applied to the simulated neutrino- induced events, and the results com- 
pared. To avoid systematic differences from data and MC, a detailed check of 
the detector simulation was performed. We used a large sample of measured 
and simulated downward going atmospheric muons, and we required that the 
distributions, used to define the selection cut, match together. The (prelimi- 
nary) results are presented in Fig. 4. It is included also the point corresponding 
to the partially contained IU events. Only 87% of IU events have been induced 



5 



1 .2 



s 1 

Q0.8 



0.6 
0.4 
0.2 




IU 



Low 



1 1 



All 



Medium 



High 



Am 2 = 2.10-"eV z 

Am 2 = 2.10" 3 eV 2 

Am 2 = 1 ,10" 2 eV 2 
maximum mixing 



27 



33 



52 

<E„> (GeV) 



Figure 4: Ratio (detected events/expected events) vs. median muon energy for the 
whole sample of 723 up throughgoing muons (All), and for the sub sample with large 
(Low, 79 events) and small (High, 122 events) deviation from a straight line. It is 
also included the ratio for the events with intermediate energy (Medium, 49 events) . 
The point IU corresponds to the 135 partially contained upgoing events. For all 
samples, the f M path length is ~ 10 4 km. The points are plotted with equal distance 
in abscissa, to guide the eye. The full line (data/MC=l) is the expectation without 
neutrino oscillations, for which a 17% total uncertainty is associated (25% for the IU 
point). The dashed lines are the expected reductions in case of — > v r oscillations 
with maximum mixing and three values of Am 2 . 



by CC interactions, the remaining by v e and NC interactions. 

Fig. 4 shows that lower energy events are more suppressed than higher 
energy ones. In a similar way, Fig. 2a shows that neutrinos with a longer 
path length are more suppressed; both effect are well explained by neutrino 
oscillation (eq. 1), with the oscillation parameter presented in Fig. 2b. 

I acknowledge all the colleagues of the MACRO collaboration, in particular 
D. Bakari and Y. Bcchcrini. 
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